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Lam inar / turbulen t transition was studied as it occurs 
in a newly constructed curved channel of rectangular cross 
section. The channel is 1.27 centimeters high by 50.8 
centimeters wide (0.5 in. x 20 in.), aspect ratio of 40 to 1, 
with a concave surface radius of curvature of 60.96 centimeters 
(24 in.). Mean velocities for air in the channel range from 
1.40 to 3.89 m/sec, which corresponds to Dean number and 
hydraulic Reynolds number ranges of 167 to 461 and 2,231 to 
6,173, respectively. 

Flow visualization from a smoke wire shows changes in the 
behavior and appearance of a smoke trace illuminated in the 
radial plane with varying Dean number and varying channel 
location. In some cases, these smoke trace patterns show 
distortion similar to what would be expected if pairs of 
counter-rotating vortices are present. A map of vortex 
stability as channel location and Dean number were varied was 
constructed from analysis of smoke trace photographs. 
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I. INTRODUCTION 



A. BACKGROUND AND LITERATURE REVIEW 

Since the early part of this century, it has been well 
known that fully developed laminar flow along a concave wall 
does not maintain a simple two dimensional (i.e. parabolic) 
velocity profile. Instead, a series of counter-rotating vortex 
pairs form as a result of varying centrifugal forces acting on 
the fluid particles. These vortices were first considered 
analytically by W.R. Dean in 1928 for flow in a narrow channel 
formed by two concentric cylinders [Ref. 1] and are similar to 
ones observed by Taylor and Gortler [Refs. 2, 3]. Figure 1 

[Ref. 4] is a schematic of Dean vortices in a curved channel. 
These vortices are known to enhance heat transfer on curved 
channel surfaces for both laminar and turbulent flows. 
Understanding such flows is expected to lead to design improve- 
ments for heat exchangers and configurations used to cool gas 
turbine blades and vanes. 

Vortex onset and development is dependent upon the magnitude 

of the curvature, as expressed by the Dean number, where 

D e = R e c y d /r Here the channel Reynolds number is defined by 

Re = Ud/v, where U is the mean streamwise velocity, d is the 
c 

channel height, r^ is the radius of curvature of the convex 
channel surface, and v is the kinematic viscosity. 
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Figure 1. Schematic of Dean vortices in a curved channel. 
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For a rectangular channel curved in the streamwise direc- 
tion, the flow at sufficiently low Dean numbers is entirely 
azimuthal. This is known as slug flow or Poiseuille flow. 

Above a critical Dean Number (De > 36.6), pure laminar flow 
becomes unstable and the secondary flow discussed above occurs. 
This is due to the fact that fluid particles in the center of 
the flow profile are travelling at higher velocities and are 
therefore subjected to higher centrifugal forces than those 
closer to the walls. The faster moving particles move outward 
toward the concave wall, forcing the fluid particles there to 
move inward toward the center of the channel where velocities 
are higher. The process repeats itself and the Dean vortices 
form. (The critical Dean number and the Dean numbers discussed 
below are the result of computer simulations using the Navier- 
Stokes equations being performed at Stanford University.) As the 
Dean number is raised further (De > 45), the vortices oscillate 
at low frequency (large wavelength) as they travel in the stream- 
wise direction. The wavelength becomes smaller up to a Dean 
number of about 75. At even higher Dean numbers (75 <06"^ 500 
to 600), unpublished simulation results (presently unconfirmed by 
experiments) show that these vortices have an oscillatory, 
rocking type of motion as they are convected. Ultimately, these 
complex secondary flow fields break down into turbulent flow. 

A brief history of earlier research into the flow and heat 
transfer characteristics of flow in a curved channel can be 
found in reference 4. The presence of Dean vortices in a 
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curved channel of aspect ratio 40 was confirmed at the Naval 
Postgraduate School by R. McKee in 1973 [Ref. 5], Modifying 
McKee’s channel and modeling it as a pair of infinite parallel 
plates with the outer wall heated and the inner wall adiabatic, 
M. Durao and J. Ballard [Refs. 6, 7] studied the effect of 

Dean vortices on heat transfer for laminar flows in straight 
and curved rectangular channels. Research continued in 1981 by 
R. Holihan, Jr. [Ref. 4] for laminar transition flows, S. 
Daughety [Ref. 8] for turbulent flows, and J. Wilson [Ref. 9] 
for transition and turbulent flows. Each of these studies 
used McKee’s apparatus, which was modified each time to improve 
the accuracy of the results. In 1985, using a new curved 
channel designed for even more precise measurements, heat 
transfer studies by G. Galyo [Ref. 10] were conducted for 
turbulent flow. 

Kelleher, Flentie, and McKee [Ref. 11] reported results of 
velocity surveys and flow visualization for the secondary flow 
in a curved channel half the size of the facility used in the 
current study. Velocity profiles were presented for Dean 
numbers of 79.2, 94.9, and 113.5. Flow visualization was 

performed using an aerosol mist (DOP) for Dean numbers of 69.3 
and 77.2. It was shown that the velocity profiles on the 
convex side of the channel mid-plane were 180 degrees out of 
phase with those on the concave side. Also shown was the fact 
that the vortices covered the full width of the channel from 
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the convex to the concave wall and that channel geometry and 
not the Reynolds number controlled vortex size. 

The numerical simulations of the Navier-S tokes 
equations for Tay lor-Gortler vortices performed by Finlay, 
Keller, and Ferziger [Ref. 12] describe changes in the vortex 
pattern which occur as the Dean number and wave number are 
varied. Also obtained are the vortex strengths and growth 
rates of the Taylor-Gort ler vortices. In addition to the 
results discussed in paragraph 3 above, this study also 
provided evidence of a small pair of secondary vortices 
developing adjacent to the convex wall for De = 80.0 and 
De = 100.0. 

B. OBJECTIVES 

The first objective of this investigation was to design 
and implement the operation of a sufficiently large curved 
rectangular duct to perform high quality flow visualization 
studies of Dean vortices. The test facility was to have 
the same aspect ratio (40 : 1) as previous investigations 

[Refs. 4 through 10], but with all dimensions scaled up by a 
factor of two. The aspect ratio is defined as the channel 
spanwise dimension divided by the channel height. 

After the channel was in operation, the objectives 
were : 
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1. to obtain information on the flow behavior at 



various Dean numbers and various streamwise and spanwise 
locations within the curved channel, 

2. to study the development of fluid path lines produced by 
a smoke wire as they are distorted and convected by the 

f lo w , and 

3. to produce stability maps for laminar-to turbulent 
transition and vortex development. 

The purposes of the results of the study are: 

1. to obtain better understanding of fluid mechanics 
phenomena which affect heat transfer in a channel where 
streamwise curvature and strong secondary flows are 
present, and 

2. to provide a basis of comparison and validation for 
computer simulations using the Navier-Stokes equations 
with appropriate boundary conditions. 

C. THESIS ORGANIZATION 

Subsequent to the introduction, this thesis describes 
the experimental facilities, detailing the specifications and 
characteristics of each component. Next, experimental 
techniques are discussed, including channel operation, flow 
visualization, photography, and measurement and calibration 
methods. Finally, experimental results and conclusions are 
presented . 
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II. EXPERIMENTAL FACILITIES 



A. CURVED CHANNEL 

A schematic of the curved channel test facility is shown 
in Figure 2. Figures 3 and 4 are front and side views of the 
actual channel. 

The lip at the entry of the inlet section is constructed of 
quarter circumference sections of six inch outside diameter 
pipe. Following this is the 25.4 cm by 50.8 cm (10 in. x 20 
in.) rectangular inlet section, consisting of a honeycomb 
followed by three screens stretched across the inlet channel 
cross section. The screens were pulled taut prior to 
tightening the flange bolts. This ensured that the screens 
were flat and perpendicular to the flow direction. The 
aluminum honeycomb is 3.81 cm (1.5 in.) deep (streamwise 
direction) and consists of hexagonal cells approximately 0.635 
cm (0.25 in.) across. The fiberglass window screen has a 
square weave with about 0.159 cm (0.0625 in.) between strands. 
These items are intended to reduce spatial non-uniformities of 
the air entering the nozzle, resulting in more stable laminar 
flow. Distance between each screen or between the honeycomb 
and the next screen is approximately 10 cm (4 in.), following 
suggestions on design rules for small low speed wind tunnels 
provided in reference 13. 
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Figure 2. Schematic of test facility. 




Figure 3. 



Test facility (front view) . 
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Figure 4. Test facility (side view). 




The air next enters a two dimensional nozzle with a 20 



to 1 contraction ratio. Designed using a fifth order 
polynomial to ensure continuous first through fourth order 
derivatives, nozzle dimensions were chosen to minimize the 
possibility of separation near the nozzle exit and to 
prevent laminar-t o-tur bulen t transition. Figure 5 shows a 
graph of the design. This nozzle induces longitudinal strain 
on the flow, stretching it to further reduce spatial non- 
uniformities which may exist. To prevent tripping the flow at 
the critical point where it enters the duct, elaborate measures 
were taken to ensure that the seam between the nozzle and the 
duct straight section was smooth. Specifically, the mating 
surface between the nozzle and duct flanges was filled with 
automotive body putty (Bondo) and sanded smooth to the touch, 
ensuring that any discontinuity (step) inside the duct was less 
than 0.001 inches. 

After the nozzle, the air passes through a 2.44 m (8 ft.) 
long straight duct. The duct inside dimensions are 1.27 cm by 
50.8 cm (0.5 in. x 20 in.). This cross section is constant 
until the air reaches the diffuser in the outlet section. 
According to reference 14, fully developed laminar flow will 
exist for hydraulic Reynolds numbers up to 6,068 (Dean numbers 
up to 454) at the end of the straight duct (i.e. entry to the 
curved channel). The corresponding mean velocity is about 3.82 
m/s. The flow is considered to be fully developed when 
normalized mean velocity profiles are self-similar . Appendix A 
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Figure 5. Graph of inlet nozzle design. 




(section 1) shows the sample calculations using the Langhaar 
equation [Ref. 14] for development length in tubes and reference 
15 for development length between flat plates. The Langhaar 
equation is shown to be the more conservative of the two. 

The fully developed laminar flow next enters the curved 
test section. Here the duct has the same rectangular cross 
section as the straight section. The concave interior duct 
surface has a radius of curvature of 60.96 cm (2 ft.) and the 
convex surface has a radius of curvature of 59.69 cm (23.5 
in.). To maintain these radii of curvature, the curved test 
section is pulled tightly into the circular support frame using 
band clamps (visible in Figure 4). Holes exist at various side 
wall locations for the installation of smoke wires. These 
wires are pulled taut in the spanwise direction in the channel. 

The air exits the curved section, goes through another 
eight foot long straight duct, and reaches the outlet section. 

At the exit of the second straight section are four screens, a 
honeycomb, a diffuser, and an outlet plenum. The screens and 
honeycomb were designed to minimize the effects of any spatial 
non-uniformities in the outlet plenum. Because of the small 
channel height, honeycomb was made of plastic straws, 0.41 cm 
(0.16 in.) in diameter and 3.18 cm (1.25 in.) long, stacked 
firmly together. The 45.7 cm (1.5 ft.) long diffuser, having a 
total angle of three degrees, is used to provide some pressure 
recovery. The design, which uses the guidance provided in 
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reference 16, minimizes the possibility of stall in the 
diffuser. The cubic outlet plenum (61 cm on each side) is 
intended to create a volume of nearly stagnant air at uniform 
pressure . 

The plenum is connected to the suction of a blower via a 
rotameter (125 cfm full scale) and a one inch standard ASME 
orifice plate shown in Figure 6. These are used for volumetric 
flow indication and calibration, respectively. Various valves 
are used to control the flow rate, and therefore the mean 
velocity and Dean number, in the curved channel. The blower 
exhaust valve is a two inch brass gate valve, while the other 
two valves shown in Figure 2 are two inch brass globe valves. 

Pipe used to connect the outlet plenum to the blower, valves, 
rotameter, and orifice is standard two inch PVC pipe. 

Both the curved and straight sections of the duct are made 
of polycarbonate (Lexan). This material was chosen because it 
is transparent and flexible, as well as being easy to join and 
machine. To minimize the potential for flow tripping, both the 
upper and lower walls of the channel sections (straight/curved/ 
straight) are made of three pieces of Lexan sheet bonded together 
in one continuous piece, so that no seams were present. Side 
wall, longitudinal, and cross-beam supports are employed, so that 
the inside dimensions of the entire channel at any strearawise 
location from the nozzle exit to the last straight section are 
approximately 0.5 + 0.005 inches by 20.0 + 0.020 inches. The 
top and bottom surfaces of the channel are secured to the side 
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Figure 6. Photograph of rotameter and orifice plate. 




walls as shown in Figure 7. Set screws in the wall segments are 
spaced about 10 cm (4 in.) apart. Wall segments are removable to 
facilitate cleaning of channel interior surfaces. Mating 
surfaces between wall segments and all flange mating surfaces 
(except for the nozzle flange) are sealed with 3-M foam tape. 
Construction of the steel support frame and all Lexan portions of 
the test facility was performed by McNeal Enterprises, Inc. of 
Santa Clara, California. 

B. BLOWER 

To provide a low pressure region at the outlet plenum, 
a Spencer VB-0370-E Vortex Blower is used. This blower is 
capable of producing about 60 inches of water vacuum at 140 
cfm and was sized by scaling flow behavior from a channel 
with half the dimensions of the present one. Figure 8 
shows the advertised blower performance curve and system 
curves. As can be seen in the figure, the blower exceeds 
the head requirements of the duct and must be throttled to 
achieve the desired flow rates. Appendix A (section 2) 
provides the sample calculations used in scaling system 
performance to determine required blower size. 

C. SMOKE WIRE 

Flow visualization is provided by smoke wires strung in the 
transverse direction at various locations in the curved 
section, as shown in Figure 9. (One wire is used at a time.) 
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Figure 7. Duct side wall design. 
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Figure 8. Blower performance and system curves. 
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Figure 9. Schematic of smoke wire, strobe, and camera placement. 



Because of the short delay between energizing the wire, smoke 
production, and arrival of that smoke a small distance down- 
stream where it is photographed, a timing control system was 
used. The control device, manufactured by Flow Visualization 
Systems of Bolingbrook, Illinois, energizes the stainless steel 
smoke wire and, after a manually selected delay, triggers the 
strobe. Voltage applied to the wire is controlled by a Calrad 
45-740 (0-130 VAC) variac. Figure 10 is a photograph of the 
control system. The delay time selected depends upon the mean 
velocity (i.e. Dean number) in the duct and was adjusted by 
trial and error. Of note is the fact that the two potentio- 
meters on the control system which control wire "on time" and 
strobe delay time do not operate independently. As shown in 
Figure 11, the sequencing circuit timing diagram, the following 
events occur during one flow visualization run: the wire is 

energized, an adjustable time delay occurs, and simulta- 
neously, the wire is turned off and the strobe is flashed. If 
the smoke wire "on time" selected is shorter than the strobe 
delay time chosen, the flash of the strobe will not occur. 

Maximum wire diameter selected is based on a wire Reynolds 

number (Re ) of 35 as discussed in reference 17. This Reynolds 
w 

number is based on wire diameter. Higher Reynolds numbers may 
result in separation and an unstable wake behind the wire. 
Sample calculations are shown in Appendix A (section 3). 

Wires penetrate the side wall segments of the test section 
through small (4.29 mm diameter) drilled holes. Pressed into 
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Figure 10. Photograph of smoke wire/strobe 
control system. 
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Figure 11. Smoke wire sequencing circuit timing diagram. 



each hole is a small ceramic bead (bead hole diameter equal to 
1.75 mm) which prevents the energized wire from burning into 
the Lexan side wall. To prevent air leakage through these 
holes, "CLAYOLA" non-hardening modeling clay is pressed flat, 
forming a patch, over each hole. Smoke wire holes are located 
at 2, 50, 112, and 140 degrees from the entry of the circular 

test section as shown in Figure 9. 

D. LIGHTING 

The strobe is a General Radio type 1531-A Strobotac, 
which is externally triggered. For the high intensity 
single flash operation used in this study, the peak light 
intensity is 7 million beam candle power. Flash duration 
is approximately 3 microseconds. The strobe, which is normally 
used to determine the speed of rotating equipment, can also 
provide more continuous illumination if flashed repeatedly at 
high frequency. The collumator lens is a 75mm (focal length), 
fl.l Rodens tock-XR-Heligon mounted directly in front of the 
strobe as shown in Figure 12. The entire assembly can be pivoted 
and clamped in place to illuminate any point on the test section. 
Since the arc produced by the strobe lamp is almost a straight 
line, the lens produces a well collumated flat sheet of light 
through the channel. 

To prevent undesired light reflections, the outer surface of 
the curved channel facing the strobe is lined with black 
construction paper (poster board). (Not shown in Figures 3 and 
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Figure 12. Photograph of strobe and collumator lens. 
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4.) Slits are cut two to ten inches downstream of each smoke 
wire to provide a path for the collimated light to illuminate a 
radial plane. To provide a basis for comparison, all slits are 
2 mm by 4.6 cm (0.08 in x 1.8 in.) and are located so that the 
center of the slit is located 4.83 cm from the center of the 
channel . 

E. CAMERA 

A Nikon F-3 with a 55mm, f2.8 lens was used to record 
the flow pattern. Kodak Tri-X (ASA 400) or Kodak Recording Film 
2475 (ASA 1000) was used. As shown in Figure 9, the camera views 
the illuminated plane of smoke particles at approximately a 
ninety degree angle to the beam of collumated light. 
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III. EXPERIMENTAL TECHNIQUE 



A. CHANNEL OPERATION 

In order to provide a variety of flow conditions for 
analysis in this study, duct flow rate (i.e. the Dean 
number) must be varied. Start-up of the channel and flow 
rate control are now described. 

To initiate flow through the duct, the blower exhaust 
and inlet bypass valves ('A' and 'B' in Figure 2, 
respectively) are verified open and the duct outlet valve ('C’) 
is checked shut. The blower is started and the duct outlet 
valve is opened. A reading of 29 percent on the rotameter can 
be achieved with all valves fully open. (Section III.B 
discusses rotameter calibration and indication.) It is 
necessary to throttle both the blower exhaust valve and the 
duct outlet valve to avoid unsteadiness (possibly due to slight 
blower stall) when operating at very low flow rates (i.e. less 
than 15% indicated flow). Failure to sufficiently throttle the 
blower exhaust valve results in severe oscillation of the 
rotameter float. Flow rates higher than 29 percent are achieved 
by throttling the blower inlet bypass valve, while valves ’A' and 
' C ' remain open. A maximum of 63 percent indicated flow through 
the duct can be achieved when the inlet bypass valve is complete- 
ly shut. Higher flow rates (up to about 80%) can be achieved if 
the orifice plate is removed. 
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B. FLOW RATE MEASUREMENTS 



Rotameter calibration was achieved by using a one inch 
diameter steel orifice, machined to ASME specifications, 
which was installed upstream of the rotameter. Orifice 
pressure taps were located at the 1-D upstream and 1/2-D 
downstream positions, where 'D* is the PVC pipe diameter 
(nominally 2 inches). Appendix A (Section 4) provides the data 
and sample calculations for rotameter calibration. References 18 
and 19 were utilized in these calculations. Data was taken at an 
ambient pressure of 30.26 inches of mercury and an air tempera- 
ture of 65°F (approx.). 

As shown in Figure 13, rotameter indication is a linear 
function of volumetric flow rate. The equations to convert 
from rotameter indication (percent) to flow rate (cfm), to 
Dean number, and to Reynolds number are shown below. 

Appendix A (section 5) provides the sample calculations for 
these results. 

Q [cfm] = 1.15 X (ROTAMETER % FLOW) + 10 

De = 9.99 X (ROTAMETER % FLOW) + 86.87 
Re c = 68.49 X (ROTAMETER % FLOW) + 595.53 

C. FLOW VISUALIZATION/PHOTOGRAPHY 

The flow visualization method is similar to that 
described in reference 20. The method described here is for 
still photography of the smoke patterns. 



37 



o 







01 
4 -> 
i — I 

P 

CO 

0) 

u 

C 



CD 

U 

P 

Cn 

*H 

iu 



38 



ROTAMETER INDICATION 



To generate an adequate amount of smoke, the smoke wires must 
be frequently wetted with smoke oil. To accomplish this, the 
smoke wires are sufficiently long at 76.2 cm (30 in.) to enable 
them to be pulled back and forth through the duct, so they can be 
wetted with smoke oil as needed. Each time the wire is used, the 
clay sealing "patches" on either wall of the duct are removed and 
the wire is slid part way out of the duct and wetted with oil 
using the fingertips or a cotton swab. Finally, the wetted wire 
is reinserted into the duct and the clay patches are reapplied. 
Variac setting for a 0.005 inch diameter wire was found to vary 
between 30 and 40 volts, depending on the channel flow rate. 

The strobe, already positioned for continuous illumination of 
the desired slit, is adjusted so that the image of the strobe arc 
is in focus on the black posterboard sheet which lines the test 
section. Focusing is accomplished by either using the adjusting 
ring on the lens or by simply moving the pivoted strobe 
reflector . 

The trial-and-error timing adjustment is accomplished by 
watching the timing control system status lights on the control 
logic box (Fig. 10). Moments after the spring-loaded timing 
sequence start switch is depressed, the smoke pattern should be 
visible in the illuminated plane. As the strobe status light 
goes out, indicating that the strobe has flashed, the smoke 
should still be visible. Adjustments to the wire and strobe 
timer potentiometers are made to accomplish this and to ensure 
that smoke intensity is at its peak when the strobe status light 
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goes out. Camera set-up and focusing is also accomplished during 
this time. 

Once the proper potentiometer adjustments are made, the 
strobe is switched to the "External Input - high intensity" 
position and photography may begin. The room is darkened as 
much as possible. (All photographs for this report were 
taken at night, since the laboratory windows could not be 
adequately shaded.) With the camera mounted securely on a 
tripod, shutter speed set to ’B', and f-stop at 2.8 (lens 
aperture fully open), the timing sequence start switch is 
depressed. Moments later, before the flash of the strobe, 
the camera shutter is opened. Immediately after the flash 
from the strobe, the shutter is closed. The focused image 
of the smoke trace appears on the film during the three 
microsecond strobe flash. 
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IV. RESULTS 



A. SMOKE TRACE VARIATION WITH DEAN NUMBER 

At a selected location in the curved channel, smoke trace 
appearance and characteristics change in a series of identifia- 
ble stages as the Dean number is increased. Results which 
follow are based on observations of smoke traces generated from 
the 0.005 inch diameter wire, and observed using continuous 
illumination over the time interval smoke was present. 
Observations are presented for different Dean numbers and 
percent rotameter flow, where rotameter full scale corresponds 
to 125 cfm at 18.3°C and atmospheric pressure. 

A summary of results in the following discussion is also 
given in Figure 15a. This figure is discussed in section IV. B. 

At the 6 degree slit location (smoke wire located 2 
degrees from the start of the curved section), the smoke trace 
appeared as a reasonably flat line for all Dean numbers from 
167 to 461 (8 to 37.5% flow by rotameter indication). For Dean 
numbers up to 407 (32% rotameter flow), smoke traces were 
particularly flat and stable. At De = 461 (approx.), the smoke 
appeared diffused across most of the illuminated plane, 
indicating that turbulence existed. Mean velocities at 
this location were the highest used, where the maximum 
corresponds to Re w equal to 31.6 for a 0.005 inch 
diameter wire at 37.5% flow (De = 461, U = 3.89 m/sec). In 
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order to avoid unstable wakes from such a wire, Re should be 

w 

less than a critical value of 35-40 [Ref. 17]. 

At the 55 degree slit location (smoke wire located 50 
degrees from the start of the curved section), the smoke trace 
appeared as a flat line at Dean numbers between 167 and 247 (8 to 
16% rotameter flow). As the Dean number was raised, the smoke 
line became blurred (De = 298 to 307). At a Dean number of 317, 
weak turbulence was observed along with a single wave peak (wave 
form) which was oscillating slightly in the transverse direction. 
At a Dean number of 327, turbulent appearance (i.e. diffused 
smoke) was observed. Due to apparent instabilities in this 
section of the duct, these effects were not completely 
reproducible each time the observations were performed: i.e., 

the same observations discussed above would occur at slightly 
different Dean numbers. During one series of observations, waves 
appeared on the smoke trace and oscillated at low frequency (3 to 
8 Hz) in the transverse direction for Dean numbers between 257 
and 307. These particular waves and wave oscillations were only 
seen during one series of observations, and were probably due to 
vibrations of the smoke wire, an effect minimized for all other 
observations . 

At the 118 degree slit location (smoke wire located 112 
degrees from the start of the curved section), observation of 
regularly spaced standing smoke pattern waves was made beginning 
with a Dean number of 167. Figure 14 shows a sequence of smoke 
pattern photographs at this location with increasing Dean 
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Figure 14a. Smoke trace photograph sequence 

(De = 182) 



II-l 2100-2330 09 FEB 1987 

9.5 % FLOW (rotameter) 

MEAN AIR VELOCITY = 1.531 m/sec 

VOLUMETRIC FLOW RATE = 0.00988 m 3 /sec 

MASS FLOW RATE = 0.01203 kg/sec 

ATMOSPHERIC PRESS. & TEMP: 29.96 in. Hg, 65.0 deg F 

SLIT LOC. = 118 deg. 

WIRE LOC. = 112 deg. 

WIRE DIA. = 0.0050 in. Re w = 19 

SLIT DISTANCE FROM START OF CURVE (X) = 1.242 m 
Re = 121,895 

De X = 182 

Re c = 1246 Re h = 2432 

KODAK TRI-X FILM ASA 400 ( f 2 . 8 , B) 
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Figure 14b. Smoke trace photograph sequence 

(De = 207) 



1 1-3 2100-2300 09 FEB 1987 

12.0 % FLOW (rotameter) 

MEAN AIR VELOCITY = 1.741 m/sec 

VOLUMETRIC FLOW RATE = 0.01123 m 3 /sec 

MASS FLOW RATE = 0.01368 kg/sec 

ATMOSPHERIC PRESS. & TEMP: 29.96 in. Hg , 65.0 deg F 

SLIT LOC. = 118 deg. 

WIRE LOC. = 112 deg. 

WIRE DIA. = 0.0050 in. Re w = 21 

SLIT DISTANCE FROM START OF, CURVE (X) = 1.242 m 
Re = 138,643 

De X = 207 

Re c = 1417 R e h = 2766 

KODAK TRI-X FILM ASA 400 (f2.8, B) 
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Figure 14c. Smoke trace photograph sequence 

(De = 227) . 



I 1-5 2100-2330 09 FEB 1987 

14.0 % FLOW (rotameter) 

MEAN AIR VELOCITY = 1.909 m/sec 

VOLUMETRIC FLOW RATE = 0.01232 m 3 /sec 
MASS FLOW RATE = 0.0150.1 kg/sec 

ATMOSPHERIC PRESS. & TEMP: 29.96 in. Hg, 65.0 deg F 

SLIT LOC. = 118 deg. 

WIRE LOC. = 112 deg. 

WIRE DIA. = 0.0050 in. Re w = 23 

SLIT DISTANCE FROM START OF CURVE (X) = 1.242 m 
Re = 152,041 

De x = 227 

Re c = 1554 Re h = 3033 

KODAK TRI-X FILM ASA 400 (F2.8, B) 
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Figure 14d. Smoke trace photograph sequence 

(De = 237) 



11-20 2100-2330 09 FEB 1987 

15.0 % FLOW (rotameter) 

MEAN AIR VELOCITY = 1.993 m/sec 
VOLUMETRIC FLOW RATE = 0.01286 m/sec 
MASS FLOW RATE = 0.01567 kg/sec 

ATMOSPHERIC PRESS. & TEMP: 29.96 in. Hg, 65.0 deg F 

SLIT LOC. = 118 deg. 

WIRE LOC. = 112 deg. 

WIRE DIA. = 0.0050 in. Re w = 24 

SLIT DISTANCE FROM START OF'CURVE (X) = 1.242 m 

Re x = 158,740 

De' = 237 

Re = 1623 Re h = 3167 

KODAK TRI-X FILM ASA 400 (f2.8, B) 
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number. (Photographs taken under other conditions and at other 
locations are provided in Appendix C.) The shape and spacing 
of the wave patterns at the 118 degree location became less 
regular as the Dean number was raised. During the three second 
time span (approx.) when the smoke trace was visible, the waves 
appeared to be motionless (De = 167 to 207). At a Dean number 
of 217, high frequency radial vibration (i.e. direction in 

Figure 1) of the smoke trace wave peaks in the illuminated 
plane began to appear. Each wave was vibrating in place at its 
own transverse location in the channel. This vibration became 
progressively more violent (higher amplitude and frequency) as 
the Dean number was raised to 247. Only the portion of each 
wave peak closest to the convex surface was observed to oscil- 
late. The approximately horizontal smoke line located between 
the convex and concave duct surfaces which connected successive 
wave peaks did not appear to be moving. (See Figure 14a.) At 
De = 252, oscillation became extremely violent and some low 
frequency transverse oscillations were evident. The turbulent 
appearance previously discussed was observed at this location 
at a Dean number of 257. 

At the 146 degree slit location (smoke wire 140 degrees 
from the start of the curved section), the smoke trace appeared 
as a series of standing waves of very regular size and spacing 
(De = 167 to 187). An observation at a Dean number of 197 
indicated that the waves had roughly the same size (amplitude), 
were spaced farther apart, and were "wobbling" in the 
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transverse (X z ) direction (i.e., the wave forms were staying 
in one transverse location, but their shape was changing in a 
wobbly fashion). This relatively low frequency wobble was 
still apparent when the smoke trace was observed at Dean 
numbers of 207 and 217, however, the wave forms were becoming 
progressively more irregularly spaced and shaped. At a Dean 
number of 227, the wobble became more violent, but the irregu- 
larly spaced and shaped waves continued to wobble in place. 

Weak turbulence (i.e., very blurred wave appearance) was 
observed at a Dean number of 237 and full turbulence appeared 
between values of 247 to 257. 

The next section on stability mapping provides graphical 
analysis of the data just described. 

B. STABILITY MAPS 

Figure 15a is a graphical presentation of the results of 
this study. As shown in the figure, the standing wave smoke 
trace was observed at the 118 degree slit location and beyond 
for the lower Dean numbers (De<210 approx.). Wave oscilla- 
tions, either a low frequency wobble or high frequency vibra- 
tion, was observed at the 55 degree slit location and beyond 
for moderate Dean values (210 < De < 310 approx.). It should be 
noted that for the 55 degree slit location, no waves were 
observed on the smoke trace until just prior to full 
turbulence. Here, the transition from weak turbulence, when 
wave oscillations were observed, to full turbulence occurred 
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Figure 15a. Flow stability map. 
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over a flow rate change from 0.02096 kg/sec to 0.02162 kg/sec, 
corresponding to a one percent rotameter flow rate change 
relative to rotameter full scale. The Dean number variation 
for this change was about 10. Turbulent smoke trace appearance 
was observed at all locations above the upper-most line in the 
figure . 

Figure 15b shows points on the stability map corresponding 
to several photographs in Appendix C. 

C. FLOW VISUALIZATION, INTERPRETATION AND ANALYSIS 

In some cases, the flow causes the smoke, which is 
initially generated along a straight line and convected down- 
stream, to deform as shown in Figure 16. One possible 
arrangement of velocity vectors in the radial plane which could 
cause such a distortion is shown in Figure 17. Figure 18 shows 
that the smoke is further distorted into shapes expected if 
streamwise vortex pairs are present. The mushroom-like shape 
shown in Figure 18 becomes more like a row of omegas as the 
smoke convects further. 

Transverse oscillations of the smoke pattern are probably 
caused by side-to-side vortex pair oscillation at a particular 
circumferential position. The causes of radial oscillations 
are less apparent. 

Analysis of the smoke patterns ran be used to make rough 
approximations of vortex strength (circula tion) and spacing. 
Figure 16 is a photograph of a smoke trace for a Dean number of 
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Figure 15b. Flow stability map indicating 
locations for several Appendix C photographs. 
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Figure 16. Smoke pattern photograph used 
in Figure 17 schematic. 



0-1 FEB 1987 

9.5 % FLOW (rotameter) 

MEAN AIR VELOCITY = 1.531 m/sec 

VOLUMETRIC FLOW RATE = 0.00988 m 3 /sec 
MASS FLOW RATE = 0.01203 kg/sec 

ATMOSPHERIC PRESS. & TEMP: NOT RECORDED 

SLIT LOC. = 119 deg. 

WIRE LOC. = 112 deg. 

WIRE DIA. = 0.0125 in. Re w = 47 

SLIT DISTANCE FROM START OF CURVE (X) = 1.256 m 

Re x = 123,238 

De = 182 

Re c = 1246 Re h « 2432 

KODAK TRI-X ASA 1,000 ( f 2 . 8 , B) 
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Figure 17. Theorized effect of streamwise vortices. 




Figure 18. Smoke pattern photograph taken farther 
downstream of Figure 17 at the same Dean number. 



11-14 2100-2330 09 FEB 1987 

9.5 % FLOW (rotameter) 

MEAN AIR VELOCITY = 1.531 m/sec 

VOLUMETRIC FLOW RATE = 0.00988 m/sec 
MASS FLOW RATE = 0.01203 kg/sec 

ATMOSPHERIC PRESS. & TEMP: 29.96 in. Hg , 65.0 deg F 

SLIT L0C. = 128 deg. 

WIRE L0C. = 112 deg. 

WIRE DIA. = 0.0050 in. Re w = 19 

SLIT DISTANCE FROM START OF CURVE (X) = 1.348 m 
Re x = 132,225 

De = 182 

Re c = 1246 Reh = 2432 

KODAK TRI-X FILM ASA 400 ( f 2 . 8 , B) 
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182 (9.5% rotameter flow). The 0.0125 inch diameter smoke wire 
was at the 112 degree location and the slit was approximately 
0.0762 meters (3 in.) downstream, which corresponds to 119.3 
degrees from the start of the curved section. Mean streamwise 
flow velocity for this Dean number is approximately 1.53 m/sec. 
(See section 5 of Appendix A for sample calculations.) Using 
this mean as an estimate of the convection velocity, the time 
between smoke generation at the wire to arrival at the illumi- 
nated plane is 0.0498 seconds. 

Assuming that the streamwise vortices (if they exist and if 
they are causing the smoke distortion) are radially centered in 
the curved channel, the maximum vortex tangential velocity can 
be calculated from an estimate of the distance from the center 
of the illuminated plane to the wave peak (i.e., wave 
amplitude). It is also assumed that the vortex is like the 
inner part of a Rankine vortex, with velocity in the radial 
plane varying linearly from the vortex center over a small 
distance. Judging from the almost linear shape of the smoke 
pattern from channel center to wave peak, this seems to be a 
reasonable assumption. Knowing that the actual channel height 
is 1.27 cm (0.5 in.) and is 4.21 cm (1.66 in.) in the photo- 
graph, measurements of smoke pattern dimensions in the 
photograph may be scaled to the correct radial lengths. In 
this case, the wave peak average amplitude is roughly 2.16mm 
from the actual channel mid-plane. This means that maximum 
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vortex tangential flow velocity must be 0.0434 m/sec, or 2.8 

percent of the mean streamwise velocity of 1.53 m/sec. 

2 

If the area, tt r , includes all vorticity of a Rankine 
vortex, then its circulation is given by T = 2irrV r , where 
V r is the circumferential component of velocity relative to the 
vortex and r is the radial distance from the vortex center. 
Using = 0.0434 m/sec and r = 2.09 mm for the smoke trace of 

Figure 16, circulation and vorticity are estimated to be 

-4 2 -1 

5.7x10 m /sec and 41.6 sec , respectively. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



Laminar/turbulent flow transition was observed in a curved 
channel of rectangular cross section. The curved test section 
consisted of concentric surfaces with radii of curvature of 
59.69 and 60.96 cm (23.5 and 24.0 in.). Channel width was 50.8 
cm and the resulting aspect ratio was 40 to 1. 

Using smoke wires at 6, 50, 112, and 140 degrees from the 
start of the curved test section, flow visualization studies 
were performed. The Dean number for these observations ranged 
from 167 to 461. Equivalent mean velocities in the duct were 
1.40 to 3.89 m/sec. 

Photographs and observations showed that under certain 
conditions and channel locations, smoke traces will appear as 
a series of standing waves at various radial planes. Smoke 
trace distortions are similar to what one would expect if pairs 
of counter-rotating vortices are present. 

Smoke pattern behavior at varying channel location and 
Dean number shows that the standing wave smoke pattern is more 
likely to form closer to the exit of the curved section at 
lower Dean numbers. As the Dean number is raised, such wave 
patterns exist farther upstream (i.e. closer to the start of 
the curvature). Transition is also expected to occur farther 
upstream, closer to the start of the curvature, as the Dean 
number increases. 
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To prove that vortices do in fact exist and that their 
presence is shown by the smoke trace, more detailed studies are 
required. The smoke trace patterns, by strict interpretation, 
show only radial distortion. Transverse flow velocities must 
also be measured to adequately show that vorticity is present. 
Radially placed smoke wires may indicate transverse flow 
characteristics, but their placement must be known relative to 
a vortex center. 

Other flow visualization and observation schemes might 
prove effective. Specifically, continuous smoke generation and 
injection of this smoke to enable following the behavior of an 
individual streamline might yield interesting results. 

The more interesting photographs taken in the course of 
this study were taken at lower flow rates and well downstream 
of the smoke wire. Since the lowest indication on the rota- 
meter was 8 percent (De = 167) and since the blower underwent 
significant throttling to achieve this flow rate, a lower 
capacity blower and smaller range rotameter would be helpful in 
further studies of transition in this test facility. 
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APPENDIX A 



SAMPLE CALCULATIONS 

1. DEVELOPMENT LENGTH FOR LAMINAR FLOW 

a. Using the Langhaar equation for tubes [Ref. 14] 

= 0.058 Re 
D c 

where 

L' = Development length for laminar flow 

D = Pipe diameter (or distance between flat plates 
(d) as discussed here) . 

Re = Reynolds number based on channel height, 

it can be shown that the maximum Reynolds number (Re ) 

max 

based on an 8 foot development length is 

L' 8 ft (12 in/ft) 

c 0.058 d 0.058 (0.5 in) 

max 

Re = 3110 

c 

max 

The maximum Dean number is 



De 



max 




3110 




0 . 5 in 



/ 23.5 



in 



De =454 
max 
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The hydraulic diameter of the channel is 



d^ = p— , where 'A' and ' P' are the 

channel area and perimeter, respectively. 

~ 4(0.5 in x 20 in) A • i_ 

D, = 7 A nn . — - — 7 0 n ex — = 0.9756 inches 
h (2x20) + (2x0.5) 

D, = 2.4 78 cm 
h 

The maximum hydraulic Reynolds number is simply 
Re,_ = Re °h = 3110 (0.9756 in.) 

h q — — 

max max d 0.5 m. 

Re, = 6068 

h 

max 

b. For development length between parallel planes, it 
was shown in Reference 15 that 
L 1 = / 

where (j) is a dimensionless parameter obtained from 
channel aspect ratio. For a channel with a 40:1 aspect 
ratio, <j) is equal to 0.0115. The maximum hydraulic 
Reynolds number is then calculated as 

_ L' _ 8 ft (12 in. /ft.) 

Ke h “ <J>D. ~ 0.0115 (0.9756 in.) 

max r h 

Re, = 8557 

h 

max 
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Similarly, the maximum Dean number and Reynolds 
number (based on channel height) using Reference 15 
results for parallel plates are 



De 



max 



Re 

c 

max 



640 

4386 



These results are about 30% higher than the more 
conservative results of the Langhaar equation for 
flow in tubes. 

2. SCALING CALCULATIONS 

To select a blower for the new test facility, which is 
twice the size of the previous facility, it is necessary to 
scale-up flow data from the smaller duct. The goal is 
having the same hydraulic Reynolds number in both channels. 
The hydraulic Reynolds number is 

Re u D h 

h . Since hydraulic diameter is 



p— and volumetric flow rate (Q) is U A, 

40 

R e h = . Ultimately, since mass flow 

rate is pQ and absolute viscosity (y) is vp, 

4m 

R e h = p^j . The perimeter (P) of the new channel is 

twice the old channel's perimeter and the cross sectional 

area (A) of the new channel is four times that of the older 

4A 

channel. Therefore, since hydraulic diameter is — , 
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D h ^ D h . Since Reynolds numbers are equal, 

new old 



similitude is achieved when mean velocities are set forth 



duct, and since the L-to-D ratio is the same for both ducts, 
similitude is also achieved when 



The equations above enable the scaling-up of flow data 
from a channel half the size of the current one. 

Figure 8 shows the old and scaled-up new curved channel 
systems curves as well as the new blower performance curve. 
Since the old channel did not have flow management devices 
(screens and honeycombs) or a rotameter, the effects of 
these items was considered small when compared to the 
suction capability of the blower and was neglected in 
generating the data for Figure 8 . 

3. REYNOLDS NUMBER BASED ON SMOKE WIRE DIAMETER 

Wire diameter used most often was 0.005 inch diameter. 
Some 0.0125 inch diameter wire was also used at lower Dean 
numbers . 



as U 




Volumetric flow rate (Q) is equal to UA. Since 



U 




new 



4 A 



old , 



Q 



new 




head loss (pressure change (AP) across the duct) is 

2 L 

proportional to U ^ , where 'L' is the length of the 



AP 



new 



equals i AP old 
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To avoid an unstable wake due to flow separation on 
the wire, maximum Reynolds number based on wire diameter 
should be less than 35 [Ref. 17] . 



(Re ) 

w 



U d 
25 max w 

v 



U (0.000127 m) 
max 

1.56x10 ^ m^/sec 



U =4.30 m/sec. (Based on Re < 35 
max w — 

for 0.005 in. diameter wire.) 

Since maximum fluid velocity for laminar flow between 

parallel flat plates is 1.5 times the mean velocity, mean 

velocity for Re^ = 35 is 2.87 m/sec. 

Reynolds and Dean numbers based on the distance between 

flat plates are now required. 

Re = Ud 
c v 

(v = 1.56x10 m^/sec) 

Re = 2,336 (Using U = 2.87 m/sec.) 
c 

The Dean number is equal to Re c , 

where d = 0.5 inches 
and r. = 23.5 inches 

l 

De = 341 (Using U = 2.87 m/sec, the maximum 

max 

velocity based on a wire diameter of 
0.005 in . ) 

(This Dean number occurs at 25.4% indicated rotameter 
flow.) For a wire which is 0.0125 inches (0.3175mm) in 
diameter, it can be shown that 
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U MAX = 1 ' 12 m/seC 
U = 1.15 m/sec 

Re =934 and 
c 

max 



De 



max 



136 



A table summarizing these values is shown below: 



SUMMARY OF CRITICAL VALUES BASED ON A WIRE 
REYNOLDS NUMBER (Re w ) OF 35 (CHANNEL INSIDE 
DIMENSIONS ARE 0 . 5in. x20 . Oin . ) : 



Wire Diameter 



U 

U 



max 

Re c 

max 



0.005 in . 
2.87 m/sec 
4.30 m/ s 
2336 



0.0125 in . 
1.15 m/sec 
1.72 m/sec 
934 



De 

max 



341 
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4. ROTAMETER CALIBRATION DATA AND SAMPLE CALCULATIONS. 



FLOW 



A one-inch ASME orifice is used, as shown above, 
is the pressure upstream of the orifice and 
is the downstream pressure. 




64 



